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CONFINEMENT IN QCD: RESULTS AND OPEN 

PROBLEMS 



Adriano Di Giacomo 



Progress is reviewed in the understanding of color confinement. 
PACS numbers: 11.15.Ha,12.38.Aw,14.80.Hv,64.60.Cn 

1. Introduction 

1.1. History 

The existence of quarks was first hypothesized by M. Gell-mann in the 
sixties 1 . The history of the idea is instructive. 

By analogy to the electromagnetic interaction it had been realized that the 
vector current of (5 decay was the Noether current of isospin conservation 
(Conserved Vector Current Hypothesis 0). 

However the electric charge Q is not a generator of the isospin symmetry 
group, since it contains an isoscalar part proportional to the hypercharge 
Y = N + S, the sum of the strangeness S and of the baryon number N. 

Q = r 3 + | (i.i) 

To put weak and electromagnetic interactions on the same footing, the sym- 
metry group had to be enlarged, to include Y among the generators, to a 
group of rank 2 containing the SU(2) group of isospin as a subgroup. Among 
the two possible candidates G2 and SU(3) the latter was found to be the 
correct choice, with hadrons assigned to the representations 1,8,10,10, the 
so called eightfold way [3]. 

All the representations of SU(3) are sum of products of the fundamental 
representations 3 and 3, so that an obvious question was about the exis- 
tence of particles in these representations, the quarks and the antiquarks, 
as fundamental constituents of hadronic matter. Their charges as predicted 
by eq(l) are fractional ±1/3, ±2/3, a clear experimental signature. 

* Presented at The Cracow School of Theoretical Physics, XLV Course. 
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An intense search for quarks was immediately started, but after 40 years no 
quark has ever been found, and only upper limits have been established for 
their production cross sections and abundance. 

It was also realized that there was a problem with Pauli principle. If the 
A 3 / 2 is made of three quarks, the state with charge 2 and spin component 
| is symmetric under exchange of the quarks since for any reasonable po- 
tential the three u quarks are in S state. A possible way out was to assign 
an extra quantum number to the quarks @], which was named color, so that 
each quark could exist in three different color states. 

After the quantization of the gauge theories it was suggested that the color 
symmetry could be an SU{2>) gauge symmetry with quarks in the funda- 
mental representation, and eight gauge bosons, the gluons mediating their 
interaction. The theory was named Quantum Chromodynamics (QCD) 5 . 
Experiments provide evidence for the existence of quarks and gluons at short 
distances, but quarks never appear at large distances as free particles. 
This phenomenon is known as Confinement of Color. 



1.2. Experiments 

The ratio R = — of the abundances of quarks and antiquarks to the 
abundance of nucleons has been investigated typically by Millikan-like ex- 
periments. No particle with fractional charge has ever been found, with an 
upper limit [Hj 

R < 1(T 27 (1.2) 

The expectation for R in the absence of Confinement can be evaluated in 
the Standard Cosmological Model |7j as follows. 

At ~ 10~ 9 seconds after Big-Bang when the temperature was T = lOGev 
and the effective quark mass m q of the same order of magnitude, quarks 
would burn to produce hadrons by the esothermic reactions 
q + q — ► hadrons 
q + q — > q + hadrons 

Putting (To = lim^o va, the burning rate is given by n q o§. The expansion 

rate in the model is equal to g]( 2 T 2 with Gn Newton gravitational constant 
ant T the temperature. The decoupling of relic quarks will occur when due 
to the burning processes the quark density will decrease to a value such that 
the burning rate is smaller than the expansion rate, or when 

n q a = G]{ 2 T 2 (1.3) 

Since the abundance of photons is ra 7 ~ T 3 , dividing both sides of eq(3) by 
T 3 gives 

^ = (1.4) 
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By use of the experimental values ^ ~ 10 9 , Gn = 10 19 /m p , and assuming 
do — "V 2 , T ~ WGev we get 

^expected = ~ = -— ~ 10 (1-5) 

Quarks have been also searched as products of particle reactions jH], again 
with no result. As an example for the inclusive cross section a q = a(p+p — > 
q(q) + X) the experimental upper limit is 

a q < 10~ 40 cm 2 (1.6) 

The expected value in the absence of confinement is &q expected — &Totai — 
10 -25 cm 2 The ratios of the upper limits to the expectations are then 

^ < ht 15 , °* < 10- 15 (1.7) 

Rexpected a 1 expected 

10 -15 is a small number. The only natural possibility is that the ratios are 
zero, i.e. that confinement is an absolute property, due to some symmetry 
of the system. 

This is similar to what happens in superconductivity, where the explanation 
for the upper limits on the resistivity is that it is exactly zero, due to the 
Higgs breaking of the conservation of electric charge, or in electrodynamics 
where the natural explanation for the upper limit to the photon mass is 
that it is exactly zero, the symmetry being gauge invariance. 
No experimental evidence exists for the confinement of the gluons. 
We shall, anyhow, define confinement as absence of colored particles in 
asymptotic states. Only color singlet particles can propagate as free parti- 
cles. 

As a working hypothesis we shall assume that some symmetry of the ground 
state is responsible for confinement. 



2. The deconfinement transition. 

A limiting temperature exists in hadron physics, known as Hagedorn 
temperature Th, due to the property of strong interactions to convert 
excess of energy into creation of particles. It was first conjectured in 1975(21 
that its existence could be the indication of a deconfining phase transition 
from hadron to a plasma of quarks and gluons. 

This transition has not yet been detected experimentally, but extensive ex- 
perimental programs and dedicated machines are being devoted to it at 
CERN SPS,at Brookhaven (RHIC),and at CERN LHC. The transition has 
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been observed in Lattice QCD. 

Both in experiments and in lattice simulations the main problem is to define 
and to detect the transition, i.e. to give an operational definition of confined 
and deconfined. In a way this problem will be the main object of my lectures. 



2.1. Finite temperature QCD 

To deal with a system of fields at non zero temperature T one has to 
compute the partition function 

Z = Tr[exp{~)] (2.1) 

with H the Hamiltonian. 

It can easily be proved that Z is equal to the Feynman Euclidean path 
integral with the time axis compactified to the interval (0, y), with periodic 
boundary conditions for boson fields, antiperiodic for fermions. 

Z = J [d<f>]e- f d3x $ £[Wf)1 * (2.2) 

A system at T = is simulated on a lattice which is in all directions bigger 
than the physical correlation length. To have a finite temperature the size 
in the time direction L t must be such that 

T = 4" ( 2 - 3 ) 
aL t 

where a = a((3, m) is the lattice spacing in physical units, which depends on 
P == ~*f and on the quark masses m. The size L s in the space directions, 
instead, must be larger than all physical scales. An asymmetric lattice is 
therefore needed L t x L z s with L s » L t . 

The dependence of the lattice spacing a((3, m) on (3 is dictated by renormal- 
ization group equations. At large enough /3's 

1 

a ^ -e 2i, o (2.4) 

with bo the coefficient of the lowest order term of the beta function, which is 
negative because of asymptotic freedom. For the temperature T of eq(2.3) 
we obtain 

A <3 

T = — el 2i, » (2.5) 
Lt 

T is an increasing exponential function of (3, i.e. a decreasing function of 
the coupling constant g 2 . This is a peculiar behavior : when the coupling 
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constant is big, and the fluctuations are large, i.e. in the disordered phase 
the temperature is small. In the ordered phase, instead, where the coupling 
constant and the fluctuations are small the temperature is large. In ordinary 
thermal systems T plays the role of the coupling constant, low temperature 
corresponds to order, high temperature to disorder. 
The key word to understand what happens is Duality. 

2.2. Duality 

Duality is a deep concept in statistical mechanics which has been ex- 
ported into field theory and string theory. 

It was first introduced in ^U] and then developed in in the frame of the 
2d Ising model which, being solvable, is a prototype system for it. 
The Ising model in 2d is defined on a simple square lattice by associating 
to each site a dichotomic field variable a = +1. The partition function is 

Z[P,a]=Y / ex P (-Y / P^j) ( 2 -6) 

ij 

The sum in the action runs on nearest neighbors and (5 = 1/T is the inverse 
temperature in units of the interaction constant. The model is exactly solv- 
able. 

A second order Curie phase transition takes place at T c = - ,„ 2 — from 

^ ^ ln(l+V2) 

an ordered ferromagnetic low temperature phase in which < a >^ to a 
disordered phase in which the magnetization vanishes. 
The model can be considered as a discretized field theory in (1+1) dimen- 
sions, and the lagrangean can be written, apart from an irrelevant constant, 
as 

L = (3 A^A.ff 

/Lt=l,2 

with AmCt = a(n + /}) — a(n). The equation of motion is A 2 a = and a 
topological conserved current exists = e fMU A u a. 

A M i^ = 

because of the antisymmetry of the tensor e^. The corresponding conserved 
charge is 



Q = = X^ e ° lAlCr ( n °' ni ) = °"( n o,+oo) - cr(n Q , -oo) 
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In a continuum version of the model, when the correlation length goes large 
compared to the lattice spacing, the value at spatial infinity being a discrete 
variable becomes a topological quantum number. Typical spacial configura- 
tions with non trivial topology are the kinks for which a is negative below 
some point n\ and positive above it. An anti-kink has opposite signs. 
It can be shown that the operator which creates a kink /i(no,roi) is a di- 
chotomic variable like a and that the partition function obeys the duality 
equation 



and the same functional form of Z on both sides of eq(2.7). 
The system admits two equivalent descriptions : 

1) a 'direct' description in terms of the fields a whose vacuum expectation 
values are the order parameters, which is convenient in the ordered phase, 
i.e. in the weak coupling regime. In this description kinks are non local 
objects with non trivial topology. 

2) a 'dual' description in which the topological excitations become local and 
the original fields non local excitations. The duality mapping eq(2.7) maps 
the weak coupling regime of the direct description into the strong coupling 
regime of the dual excitations and viceversa. The dual description is con- 
venient in the strong coupling regime of the direct description. 

The 2-d ising model is self-dual, being the form of the dual partition func- 
tion the same as that of the direct description, but this is not a general fact. 
Other examples of duality are : the duality angles-vortices in the 3-d X-Y 
model ^21) the duality magnetization - Weiss domains in the 3d Heisen- 
berg model ^3], the duality ^-monopoles in compact U(l) gauge theory 
dUUSlCni) the duality fields-monopoles in N=l SUSY SU(2) gauge theory, 
and many examples in string theory ^SJ. 

The idea is then to look for dual, topologically non trivial excitations in 
QCD, which we shall generically denote by fi, which are ordered in the 
confining phase < fj, >^ 0, thus defining the dual symmetry. 



The same problem as in experiments exists for Lattice simulations : how 
to define and detect the confined and the deconfined phase. 
In pure gauge theory (no quarks, quenched ) the Polyakov criterion is used, 
which consists in measuring the qq potential at large distances. If it grows 
linearly with distance 



Z[f3,a] 



Z[f3*,fj] 



(2.7) 



with 




2.3. The deconfinement transition on the Lattice 
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V{R) R ^oo ~ aR (2.8) 
there is confinement. If it goes to a constant 

V(R) R -> 0O ~C+ < ^ (2.9) 

the phase is deconfined. The potential is measured through the correlator 
of Polyakov lines. A Polyakov line is defined as the parallel transport along 
the time axis across the Lattice. 

L(x) = Pexp( igA (x,t)dt) (2.10) 
Jo 

In terms of the correlator of two Polyakov lines 

G(x-y) =< L(x)L(y) > 

the static potential V(x — y) acting between a quark and an antiquark is 
given by 

V(x-y) = -Tln{G{x-y)) (2.11) 
At large distances, by cluster property, 

< L(x)L(y) >| 5 f_^- oo| w | < L > | 2 + Kexp(- ^ X ~^ ) (2.12) 

If | < L > | ^ then V(R) — ► constant as R — > oo and there is no 
confinement. If, instead, | < L > \ = then, at large R, V(R) ~ crR and 
there is confinement. 

| < L > | is an order parameter for confinement, Zn, the centre of the gauge 
group, being the relevant symmetry. 
Indeed it can be shown that 

\<L>\= exp(-^) (2.13) 

with F q the chemical potential of a quark. In the confined phase F q diverges 
and | < L > \ —> 0. 

There is a problem in the continuum limit since F q diverges also in the 
deconfined phase due to the self-energy of the quark, and a renormalization 
is needed [TU]. 

A transition is observed on the Lattice at a temperature T c from a low 
temperature phase where | < L >= 0| (confinement) to a high temperature 
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phase where | < L >^ 0| (deconfinement). 

For gauge group SU(2) = .50 and the transition is second order in the 
universality class of the 3d ising model. 

For gauge group SU(3) ^ = .630(5) and the transition is weak first order. 
With the usual convention 2ira = lGev this gives T c ~ 270Mev. The order 
of the transition is determined by use of finite size scaling techniques, which 
are nothing but renormalization group equations [See e.g. |20|]. 
The density of free energy F by dimensional arguments depends on the 
spacial size L s of the system in the form 

/ = F(L s )Lj = *(° f ) (2.14) 

where a is the lattice spacing and £ is the correlation length. 

In the vicinity of T c £ goes large with respect to a, so that | « 0. 

Since £ diverges as r = (1 — — > as 

foci- - " (2.15) 
the variable 4- can be traded with the variable tL's , and 

f = <KtL») 



For the specific heat Cy = —y-^j and for the susceptibility x<L> 
j d?x < L(x)L(0) > the resulting scaling laws are 

C v -C = Lf^c(TLl) (2.16) 
X< l> = lU <l> (tlI) (2.17) 

From the measured behavior with L s of these quantities the critical indexes 
a, 7, v can be determined, which identify the universality class of the tran- 
sition. 

For 3d ising a = .11,7 = 1-43, v = .63. 

For a weak first order a = 1,7 = 1, v = 1/3. 

In the presence of quarks Z3 is not a symmetry any more and < L > is not 
an order parameter. The string breaks also in the confined phase and its 
energy is converted into pions. How to define confined and deconfined? 
The phase diagram for the case Nf = 2 with two quarks of equal mass m 
is shown in Fig^ A line exists across which < L >, < ipip >,< E > all 
experience a rapid change, so that their susceptibilities have a peak. All 
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Fig. 1. The Phase Diagram of JV> = 2 QCD 



these peaks happen to coincide within errors. Conventionally the phase 
below the line is called confined, the one above it deconfined. An order 
parameter is needed, which must exist if, as we have argued, the transition 
is order-disorder. The dual excitations have to be identified. 

3. The dual excitations of QCD 

The general idea is that the low temperature phase of QCD (strong cou- 
pling) can be described, in a dual language, in terms of topologically non 
trivial excitations which are non local in terms of gluons and quarks, but 
are local fields in a dual language, and weakly coupled [2~T] [2*2*] [T7j . 
There exist two main proposals for these excitations, both due to 'tHooft. 

1) Vortices PT] 

2) Monopoles E2,|23| , EH- 

In the vicinity of the deconfining transition (T < T c ) the free energy 
density should depend on the dual fields in a form dictated by symmetry 
and scale invariance. The deconfining transition is a change of symmetry : 
the disorder parameter < $duai for T < T c , < $duai >= for 

T > T c . 

Two main approaches have been developed in the literature: 

a) Expose in the lattice configurations the dual excitations and show that 

by removing them confinement gets lost. (Vortex dominance, abelian dom- 
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inance, monopole dominance) 

b) Study the symmetry and the change of symmetry across T c . 



3.1. Vortices 

Vortices are one dimensional defects associated to closed lines C, V(C). 
If W(C) is the Wilson loop, i.e. the parallel transport along the line C 
then 

V{C)W{C) = W{C')V{C)e 1 ^^ (3.1) 

ncc' is t ne linking number of the two curves C, C , which is well defined in 
3d. 

In (2 + l)d, < V(C) >^ means spontaneous breaking of a symmetry, the 
conservation of the number of vortices minus the number of anti- vortices, 
< V(C) >= means super-selection of that number, and V(C) can be an 
order parameter for confinement. In (3+l)d this statements have no special 
meaning. 

In any case, as a consequence of Eq(3.1) whenever V(C) obeys the area law 
W(C) obeys the perimeter law, and viceversa whenever W(C) obeys the 
area law V(C) obeys the perimeter law. 

The 'tHooft loop, defined as the expectation value of a vortex going straight 
across the lattice, or the dual of the Polyakov line, is non zero in the confined 
phase, zero in the deconfined phase :25|. The corresponding symmetry is Z3, 
which, however, does not survive the introduction of dynamical quarks. 



3.2. Monopoles 

Monopoles exist as solitons in Higgs gauge theories with the Higgs in 
the adjoint representation |26| |27| . They are stable for topological reasons. 
If the gauge group is SU(2) they are hedgehog-like configurations for the 
Higgs field (/>, with 4> l (f) cx r l , and are characterized by a zero of <p corre- 
sponding to the position of the monopole. These configurations are called 
monopoles because of the non trivial topology of the mapping of the sphere 
at spacial infinity S2 on the sphere of the possible values of < cp >. 
Physically this Ccin be understood, in terms of the 'tHooft tensor, F^ u . 

F„ u = ^G^ - -4>(D U 4> A D v 4>) (3.2) 
9 

where g is the gauge coupling constant, D^4> the covariant derivative of (j), 
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and^. 

F^ u is gauge invariant by construction. Moreover the bilinear terms in A^A U 
cancel between the two terms of eq(3,2) and 

In the unitary gauge <j) = (0, 0, 1), the last term vanishes and 

Ffiv = 9/j.A^, — dyA"^ 

an abelian field. 

F^ v obeys Bianchi identities 

d„F; v = (3.3) 

with F* u = ^e^ppFpo- the dual tensor. The identity can be violated at the 
location of singularities, where a non zero magnetic current exists 

In any case due to the antisymmetry of F* u 

d v J v = (3.4) 
For the monopole solution j^f)] 

FfiQ = 0,(E = 0) 



i i n 

2 e ijk F jk = 2^ + Dirac - string 



A Dirac monopole. The string is produced by the singularity of the 
transformation to the unitary gauge at the zero of <j>. The transformation 
to the unitary gauge is called Abelian Projection. 

For SU(N) gauge group one can inquire about the existence of monopole 
solitons and what the analog of (j) is[2Hl- ^ we denote by $ = Y, a & a T a the 
Higgs field, by A^ = S a ^^T a the gauge field and by = d^A v - d v A^ + 
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ig[A^,A u ] the field strength tensor, with T a the generators of the gauge 
group in the fundamental representation, normalized as Tr(T a T b ) = 5 ab ,we 
can define the generalized 'tHooft tensor as 

F/J.U = Tr(<Z>G^ - ^[D^, A,*]) (3.5) 

The necessary and sufficient condition to have abelian projection, i.e. can- 
cellation of bilinear terms in A^Ay is that 

* = <f>*,<f>" = U(x)*$% ag U(x) 



with U(x) an arbitrary gauge transformation and 

,- < a a a(N-a)times (N-a) (N - a) (N-a). 1 n , A _ _ 

(3.6) 

The residual symmetry is SU(a) <g> SU(N - a) ® U(l). 
For each <E> a one has 

F« v = Tr(<S> a G^--<S> a [D ll <S> a ,D u <S> a ]) = d^Tr(<S> a A u )-d u Tr($ a A f ,)--Tr(<5> a [d^ a ,d u <5> a ]) 

5 (3.7) 

Transforming to the unitary gauge where $ a = <&% ag gives 

F« v = d^Tr(^ a diag A u ) - dvTr^AJ (3.8) 

Expanding the diagonal part of A^ as a sum of simple roots of the algebra 
of the group, a a , which obey the orthogonality relations Tr(a a <& b diag ) = 5 ab , 

A dia 9 = £ a AX,one g ets 

Kv = d »K - duAl (3.9) 
which is an abelian field. The simple roots have the form 

a a = dia0(O,O,....O,l,-l,O,...O) 



with the 1 at the a-th entry. A monopole soliton solution exists for each 
value of a in the SU{2) subspace spanned by the elements +1 and —1. For 
the Higgs field one has 

<&(*) = U(x)^(x) diag U(x) 
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where <&{x)diag is defined with eigenvalues in decreasing order. Expanding 
${x) diag in the complete basis $ diag , 

${x) diag = Z a c a (x)$ a diag 

one gets 

= £ a c a (x)$ a (x) (3.10) 

The transformation U(x) is singular at the sites where some c a {x) vanishes, 
i.e. wherever two subsequent eigenvalues of $ coincide: these points are the 
locations of the monopoles. The field strength F" can be defined also in 
the absence of a Higgs field in the lagrangean simply as 

¥%, = Tr($ a G^ - -^ a [D^ a , D u <S> a \) (3.11) 
^(x) = U\x)^ diag U{x) (3.12) 

U(x) an arbitrary gauge transformation which can have non trivial topology 
or singular points. 

F£ v depends on the choice of U(x). It obeys the Bianchi identities 
Eq.(3.3), apart from singularities where the magnetic current can be non 
zero. 

In any case the magnetic current will obey the conservation law eq(3.4). 
The theory has (N-l) topological symmetries built in, corresponding to the 
conservation of magnetic charges. 

If these symmetries are realized a la Wigner the Hilbert space will be supers- 
elected. If they are Higgs broken the system will be a dual superconductor. 
Our working hypothesis will be that the dual symmetry of QCD is the 
conservation of (N-l) magnetic charges. The change of symmetry at T c is a 
transition from Higgs-broken to superselected. Dual excitations carry mag- 
netic charge. 

Our program will then be to construct magnetically charged operators \i a 
and study their vacuum expectation values < \x a >. 

< n a >^ means dual superconductivity. 

< fi a >= means normal vacuum. 

This should hold both in quenched theory and with dynamical quarks, in 
agreement with the ideas of iV c — > oo limit of QCD. 

3. 3. Construction of < /i a > 
The basic idea is simply that 



e ipa \x >= \x + a > 



(3.13) 
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if x is a position variable and p its conjugate momentum. 
Specifically 

M a (x,i) = exp(i J d 3 yTr[$ a (y,t)E(y,t)]b ± (x-y)) (3.14) 

where E is the electric field operator, and b(x — y)± is the vector potential 
produced by a static monopole sitting at x in y. 

Vb ± = 0, 
-* j* 2-k f 

V A b± = t + Dirac — string 

9 r 6 

**{x) = rf(x)*% ag U(x) 

with U(x) a generic gauge transformation. 
u a is gauge invariant. In the gauge 

$ a = $5 

a a {x,t) =exp(i J Tr{$% ag E(y,t)]b ± {x-y)d 3 y = exp(i J E a ± (y,t)b ± (x-y)d 3 y 

(3.15) 

Ej_(y,t) is the conjugate momentum to Aj_(y,t) so that 

fi a (x,t)\A\(y,t) >= \A a ± (y,t) + b ± (x-y) > (3.16) 

A Dirac monopole has been added to the abelian projected configuration. 

There are (N-l) species of monopoles, corresponding to a=l, N-l. 

a a creates a singularity (monopole) in a selected gauge and in all the gauges 
obtained from it by a transformation which is continuous in a neighborhood 
of the singularity. The number of monopoles per fm 3 is finite as illustrated 
in Fig. 2 where an histogram is displayed of the distribution of the difference 
between two eigenvalues of a plaquette operator for different values of the 
lattice spacing 29. . Therefore creating a monopole in an abelian projection 
implies that a monopole is also created in any other abelian projection, apart 
from a set of zero measure. The statements < a a >= and < u a ^ > 
are independent of the abelian projection, so that the statement that QCD 
vacuum is or is not a dual superconductor are absolute, projection indepen- 
dent statements. 
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— Lattice 32 
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= 6.6 
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Lattice 24 
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= 6.4 


4 

— Lattice 16 


p= 


= 6.0 



Fig. 2. Distribution of the differences of the phases of the eigenvalues of the 
Polyakov line, for three lattices with the same physical volume and different lattice 
spacing. A monopole on any site would correspond to a non zero value at zero 
angle. 



3.4- Measuring < p a > 

By construction 

< ^ >= ^ (3-17) 

where Z is the partition function of the theory, and Z a the one modified by 
the insertion of the monopoles. Eq(3.17) implies that at (3 = < p a >= 1. 
Taking advantage of that it is convenient, instead of measuring < p a > 
directly, to measure its susceptibility p a = -§pln < fi a > which is much less 
noisy and will prove more suitable for our purposes. From Eq(3.17) one 
immediately gets 

p a =< S > s - < S a > s « (3.18) 

One also has 

< fi a >= exp[ J dp'p a (P')] (3.19) 

It follows from Eq(3.19) that, in the infinite volume limit : 

(i) < p a for T, T c iff p a tends to a finite limit. 

(ii) < p a >= for T > T c iff p a — > oo The property (ii) is much easier to 
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check on p than by a direct measurement of < p > which can only give 
limits, due to statistical errors. 

In the critical region T ~ T c a strong negative peak is expected due to a 
rapid decrease of < p a >, and scaling laws corresponding to the fact that 
the correlation length goes large with respect to the lattice spacing. The 
renormalization group equations read 30 

< p a >= L1$ a (rLS ,mL* h ) (3.20) 
i 

Sending L s — > oo keeping tL% fixed gives 30 

-K 1 

< p a >^mn cj) a (TLf) 

or 

f/LS = f(rLX) 

a scaling law from which the critical index v can be determined. 
The prototype theory is compact U(l) in Ad, where everything is understood 
analytically at the level of theorems ^1] ^3] |31| . There is a phase transition 
at P c ~ 1.01 which is first order, from a confined phase to a deconfined phase, 
and < p > is non zero below f3 c and zero above (5 C . 

Moreover p is proved to be a gauge invariant charged operator of the Dirac 
type. A numerical determination provides a check of the approach. The 
result is shown in Fig.sEGEl 

A strong negative peak signals the transition. At low fi' s p is size in- 
dependent, at large (3' s it is proportional to L s with a negative coefficient, 
implying that < p > is strictly zero in the thermodynamic limit. 

Finite size scaling agrees with a first order transition. 

For quenched SU(2) theory the deconfining transition is detected in a 
similar way at the right value of (3 and the critical index v is that of the 3d 
ising model |3*2]. 

Quenched SU{3) also shows a first order transition at the right tem- 
perature |3"3] . 

The numerical check of the independence on the choice of the abelian 
projection is contained in |34j . 

The case of Nf = 2 QCD can be approached in the same way. The 
results are displayed in the Fig.'jol 17151 1§1 The finite size scaling is that of 
a first order transition, and definitely excludes a second order transition in 
the universality class of 0(4), 0(2) model[2*9]. 

The implications of this fact, together with a finite size scaling analy- 
sis of other quantities, like the specific heat, the chiral condensate and its 
susceptibility will be the object of the next section. 
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Fig. 4. Size dependence of p below (3 C 



4. Nf = 2 QCD 

QCD with two flavors of light quarks is a good approximation to nature, 
and also a specially instructive system from the theoretical point of view. 
For the sake of simplicity we shall consider two quarks of equal mass m. 
The phase diagram is shown in Fig^ For m > 2Gev the system is quenched 
to all effects, the phase transition is first order and < L > is a good order pa- 
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rameter, Z% the relevant symmetry. At m ~ a phase transition takes place 
from the spontaneously broken phase to a symmetric phase, and < i()tp > is 
the order parameter. In the intermediate region of m's chiral symmetry is 
broken by the mass, Z3 is broken by the coupling to quarks and apparently 
there is no order parameter. Also The Ua(X) symmetry which is broken 
by the anomaly is expected to be restored about at the same temperature 
as the chiral symmetry. Three transitions, deconfinement, chiral, Ua[X)' 
are they independent? Of course a definition of deconfinement is needed to 
answer this question. 



4-1. The Chiral Transition 

If one assumes, following reference |35| that low mass scalars and pseu- 
doscalars are the relevant degrees of freedom, the order parameters are 

* : $ij =< ^(1+75)^ >,i,j = 1..JV) (4.1) 

Under the symmetry group SU(Nf) ® SU(Nf) <S> Ua(X) 

$ -► e ia U L $U R 
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. The most general effective Lagrangean (density of free energy) invariant 
under the symmetry group is 

i 2 2 2 

L eff = -Tr(5 M $t^$)_!^ Tr ($t$)_^_ 5l [ Tr ($t$)]2_l_ 52 [ rr ,^t$)2] +c [ de ^ +de ^t] 

(4.2) 

Terms with higher dimension have been neglected since they become irrel- 
evant at the critical point. 

Infrared stable fixed points indicate second order phase transitions. A (4— e) 
extrapolation to 3d is intended. The last term in Eq(4.2) is the Wess-Zumino 
term describing the anomaly: it is invariant under SU(Nf) <8> SU(Nf) but 
not under Ua(1), and has dimension Nf. For Nf > 3 it is irrelevant and no 
IR stable fixed point exists, so that the transition is weak first order. 
For Nf = 2 instead the Wess-Zumino term has dimension 2 so that its 
square and its product with the mass term are also relevant. If c = at the 
fixed point the symmetry is 0(4) ® 0{2) and no IR fixed point exists, so 
that the transition is 1st order. Physically this happens if the mass of the 
r)t, m v i, vanishes at T c . 

If c ^ 0, or if m^i is non zero at T c , the symmetry is 0(4) and the transition 
can be second order. If this is the case the transition is a crossover around 
the critical point [see Fig. a tricritical point is expected at non zero 
chemical potential jHB] which could be observed in heavy ion collisions. No 
evidence of it has emerged from experiments to date. 

If, instead, the transition is first order it will also be such in the vicinity of 
the chiral point and possibly all along the transition line, and no tricritical 
point exists. 

This issue is fundamental to understand confinement: a first order phase 
transition is a real transition and can correspond to a change of symmetry 
and to the existence of an order parameter. 

A crossover means that one can go continuously from the confined region 
to the deconfined one and that confinement is not an absolute property of 
the QCD vacuum. 

4-2. Thermodynamics 

The order of the transition can be determined by a finite size scaling 
analysis |37| |38l of lattice simulations. 

Let r = (1 — tjt) be the reduced temperature. As r — > the correlation 
length of the order parameter, £, diverges as 

so that the ratio of the lattice spacing a to £ is negligible and there is scaling. 
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If L s is the spacial extension of the lattice, the scaling laws read 

C v - C « Lf^>c(rLl,mLf) (4.3) 

and 

X nLh x (TLg,miP) (4.4) 

Here Cy = f^y, and 

X = J d 3 x< $(s)*(d) > co?m 
is the susceptibility of the order parameter 

The critical indexes a, /3, 7, f are the anomalous dimensions of the operators 
and identify the order and the universality class of the transition. Eq's(4.3) 
and (4.4) are nothing but the renormalization group equations. The sub- 
traction needed for Cy corresponds to an additive renormalization 38 . 
Notice that the scaling law of the specific heat is unambiguous whilst that 
for x only holds if $ is the order paramenter : the equality of the index v 
for the two scalings can a legitimation of the order parameter. 
The scaling laws Eq.'s(4.3) and (4.4) involve two scales, a fact which makes 
the analysis complicated with respect to the simpler case of quenched QCD. 
To simplify the problem one can study the dependence on one scale by keep- 
ing the other one fixed |3U| . One possibility is to vary m and L s keeping the 

quantity mLg which appears in the scaling laws fixed. The scaling equation 
(4.3) becomes then 

C v - C ~ L£$ c (tlI ,mV s h = M) (4.5) 
so that the peak scales as 

{C V - Co)peak oc Lf (4.6) 

This allows a determination of ^. The critical index yh is the same within 
errors for 0(4) and 0(2) universality classes, so that the same simulations 
can be used to check both universality classes; moreover the index a is 
negative for both, implying that the peak should decrease with increasing 
volume, (see Table 1) FierlTUlshows a test of Eq(4.5),FigH3a test of Eq(4.6). 
0(4) and 0(2) universality classes are excluded with a high confidence level 
(x 2 /dof ~ 20) : the peaks increase rapidly with the volume instead of 
decreasing. 

For the details of the determination of the subtraction Co see [30] . 
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Table 1. Critical exponents. 
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Fig. 10. Test of scaling Eq.(4.5), with M = 74.7 (left) and M = 149.4(right) : the 
curves should coincide. 



A similar result is obtained for the susceptibility of < ipip > which is 
believed to be a good order parameter near T c (Fig^Jand Fig. [T^]l 

X ~ $ x (TLj,mL y s h = M) (4.7) 

and 

Xpeak oc Lj (4.8) 

For this test x 2 /dof w 10. 

As a result we can state that the transition is neither in the universality 
class of 0(4) nor in that of 0(2). Another possibility is to look at the large 




Fig. 11. Scaling Eq.(4.7) of the chiral susceptibility, with M = 74.7 (left) and 
M = 149.4 (right) : the curves should coincide. 
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Fig. 12. Test of the scaling Eq.'s (4.6) and (4.8) 



volume limit keeping the first variable fixed: tLs is related to the ratio -g- 
of the correlation length to the spacial size of the lattice, while the other 
variable mL y s h is related to the ratio of the pion Compton wave length ^- 
to L s . As L s goes much larger than ^- a finite limit is reached and [30] 

C v -C Btm^fc(TLz) (4.9) 

and 

7 1 

Xttm v *h$ x (rLZ) (4.10) 

The result of this analysis is shown in Fig ll3b nd Fi dl4l No scaling is 
observed assuming second order transition with 0(4) or 0(2), but a good 
scaling for first order. 

A similar result is obtained for the scaling Eq.(4.10) Fig |15l 
Finally one can investigate the so called magnetic equation of state: 

i -i 

< ipip >= ms f(Tm" y h) (4.11) 

For 0(4) 5 = 4.85, for first order 5 = oo Again no scaling is observed 
assuming 0(4) or 0(2) second order transition Figll6| and good scaling for 
first order, Fi dTTl 

The issue is fundamental and deserves further attention. 



TEMPLATE PRINTED ON FEBRUARY 1, 2008 





100 



Fig. 14. Scaling Eq.(4.9) for first order. 



template printed on February 1, 2008 



27 




Fig. 15. Testing Eq.(4.10) for 0(4) and first order. 
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Fig. 16. Scaling of the magnetic equation of state (4.11) assuming 0(4) 



5. Concluding remarks 

We have discussed the experimental evidence for confinement and how it 
naturally implies that there exists a dual symmetry in QCD whose breaking 
is responsible for confinement. We have presented the two most accredited 
candidates for dual topological excitations, vortices and monopoles. 
We have then shown how the working hypothesis that monopoles confine 
via dual superconductivity of the vacuum can be tested by numerical sim- 
ulations on the lattice, through an order parameter which is the vev of 
operators carrying magnetic charge. The numerical tests strongly support 
the validity of this idea, which can be put in a consistent form and made 
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Fig. 17. Scaling of the magnetic equation of state (4.11) assuming first order 



independent on the choice of the abelian projection. This holds both for 
quenched QCD and in the presence of dynamical quarks. 
A prerequisite is that deconfinement is a true order-disorder phase transi- 
tion, and not a crossover, which would allow a continuous path from confined 
to deconfined phase. We have thus discussed a test with Nj = 2 QCD where 
an unsolved dilemma exists between the existence of a crossover and a first 
order phase transition. We definitely exclude a second order chiral transi- 
tion which would imply a crossover at non zero quark mass, whilst we find 
evidence for a first order transition. The issue is fundamental and deserves 
further studies. 

From what we have seen we can conclude that the dual excitations of QCD 
are magnetically charged, or that dual superconductivity of the vacuum can 
be the mechanism of confinement. 
However we are not yet able to identify them. 
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